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Observations of gravitational waves from inspiralling compact binaries using laser-interferometric detectors can provide accurate 
measures of parameters of the source. They can also constrain alternative gravitation theories. We analyse inspiralling compact 
binaries in the context of the scalar-tensor theory of Jordan, Fierz, Brans and Dicke, focussing on the effect on the inspiral 
of energy lost to dipole gravitational radiation, whose source is the gravitational self-binding energy of the inspiralling bodies. 
Using a matched-filter analysis we obtain a bound on the coupling constant ubd of Brans-Dicke theory. For a neutron- 
star/black-hole binary, we find that the bound could exceed the current bound of ojbd > 500 from solar-system experiments, 
for sufficiently low-mass systems. For a 0.7Mq neutron star and a 3M© black hole we find that a bound uibd ~ 2000 is 
achievable. The bound decreases with increasing black-hole mass. For binaries consisting of two neutron stars, the bound is 
less than 500 unless the stars' masses differ by more than about 0.5Mq. For two black holes, the behavior of the inspiralling 
binary is observationally indistinguishable from its behavior in general relativity. These bounds assume reasonable neutron-star 
equations of state and a detector signal-to-noise ratio of 10. 



PACS numbers: 04.30.-w,04.80.Cc 
Preprint No. WUGRAV 94-6 

I. INTRODUCTION AND SUMMARY 

The regular detection of gravitational radiation from 
astrophysical sources by large-scale laser-interferometer 
systems as the US LIGO or the French-Italian VIRGO 
projects will usher in a new era of gravitational-wave as- 
tronomy jjj . One of the most promising sources for de- 
tection by laser-interferometric systems is the inspiralling 
compact binary, a binary system of neutron stars or black 
holes whose orbit is decaying toward a final coalescence 
under the dissipativc influence of gravitational radiation 
reaction. For much of the late-time evolution of such 
systems, the gravitational waveform signal is accurately 
calculable 0), given by a "chirp" signal, increasing in 
amplitude, and sweeping in frequency through the de- 
tectors' typical sensitive bandwidth between 10 Hz and 
1000 Hz. Astrophysical estimates of the rate of such in- 
spiral events are promising: for the advanced version of 
LIGO, capable of detecting the inspiral waveform to dis- 
tances of hundreds of Mpc, the estimated rate is 3 per 
year, and could be as large as 100 per year ||. 

In addition to simple detection of the waves, it will be 
possible to determine important parameters of the inspi- 
ralling systems, such as the masses and spins of the bod- 
ies Q . This is made possible by the technique of matched 
filtering of theoretical waveform templates, which depend 
on the system parameters, against the broad-band de- 
tector output H|. The method exploits the fact that, 
depending on the source, between 500 and 16,000 cycles 
of the waves may be observable in the sensitive band- 
width, and so the matching of a template to the sig- 
nal will be extremely sensitive to the evolution of the 
gravitational-wave frequency with time. That evolution 
depends, of course, on gravitational radiation reaction, 
which depends on the parameters of the system. Very 



accurate determinations of the masses of the components 
should be possible, while less accurate estimates of spins 
and other parameters may be feasible f|[jj- 

It is said that the first detection of gravitational radi- 
ation will also constitute a verification of general relativ- 
ity, since that is the basic theory used in all calculations 
of gravitational radiation from such systems. It would 
be useful, however, to quantify that statement, by as- 
sessing how accurately such observations could actually 
constrain or bound alternative theories of gravity. This is 
not a straighforward question to answer with generality. 
Whereas in the slow-motion, weak-field, non-radiative 
limit appropriate to solar-system dynamics, most alter- 
native metric theories of gravity can be encompassed 
by one simple framework, known as the parametrized 
post-Newtonian (PPN) formalism (see Q for review and 
references), no correspondingly simple framework exists 
for describing radiative systems, or systems containing 
strong-internal-field, compact objects, such as neutron 
stars or black holes. On the other hand, for one simple, 
but popular class of alternatives, the scalar-tensor theory 
of Fierz, Jordan, Brans and Dicke ipyTof, and some of its 
generalizations, the full details have been worked out. In 
this paper, we shall explore the extent to which obser- 
vations of gravitational waves from inspiralling compact 
binaries could usefully constrain scalar-tensor gravity. 

For simplicity, we focus on the version of scalar-tensor 
gravity known for short as the Brans-Dicke (BD) the- 
ory. That theory augments general relativity (GR) by 
the addition of a scalar gravitational field that couples 
universally to matter (hence the theory, like GR, is a met- 
ric theory, satisfying all fundamental equivalence princi- 
ple tests |ll]]) and determines the gravitational coupling 
strength G via G oc (f> . The relative importance of the 
scalar field is parametrized by a coupling constant Wbd 
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(in generalized scalar-tensor theories, wbd can itself be 
a function of the scalar field) . Roughly speaking, in the 
limit of large wbd, the relative difference between effects 
in GR and effects in BD is 0(1/wbd)- As wbd — ► oo, BD 
tends smoothly toward GR. The best current empirical 
bound is cjbd > 500, from solar-system measurements of 
the Shapiro time delay and the deflection of radio waves 
by the Sun f§. 

For systems involving gravitational radiation and com- 
pact objects, BD introduces three effects [ ^3| , ^4|Jl^| : 

(i) Modifications to the effective masses of the bod- 
ies. These modifications depend on the internal struc- 
ture of the bodies, as parametrized by "sensitivities" sa, 
which roughly measure the gravitational binding energy 
per unit mass. These effects violate the Strong Equiv- 
alence Principle in that the motion of such bodies 
now depends on their internal structure (apart from tidal 
interactions). For neutron stars, s ~ 0.1 — 0.2, and for 
black holes s = 0.5. 

(ii) Modifications of quadrupole gravitational radia- 
tion. BD predicts monopole as well as quadrupole grav- 
itational radiation, whose combined effect is to modify 
the effective GR quadrupole formula for two-body energy 
loss, 
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by corrections of C(1/wbd)- Here, p, and m are the re- 
duced and total mass respectively; r, v and r are the 
orbital separation, velocity and radial velocity, and the 
units are such that G — c = 1 . 

(iii) Dipole gravitational radiation. The center of grav- 
itational binding energy need not be coincident with the 
fixed center of inertial mass, if the two bodies are differ- 
ent, and in BD, the resulting varying dipole moment is a 
source of scalar radiation. Because it is a dipole rather 
than a quadrupole effect, the dipole contribution to the 
energy loss has two fewer time derivatives, and thus is 
0(v~ 2 ) larger than the quadrupole contribution. It also 
depends on the difference in sensitivities, <S = Si — s%, 
between the two bodies. Specifically, 
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We work here to first order in 1 / wbd • 

The most important consequence of dipole gravita- 
tional radiation is that it modifies the evolution of the or- 
bital radius and thence the gravitational-wave frequency 
/, because /// = -(3/2)r/r = -(3/2)E/\E\. In the 
matched filtering method, any difference between the fre- 
quency evolution of the theoretical template and that of 
the actual signal will ultimately cause the two to go out 
of phase, and the signal-to-noise ratio will drop. A rough 
measure of the accuracy of the template, then, can be 
obtained by determining how much of a change in the 
template is sufficient to cause a change of ~ ir radians in 



the total accumulated gravitational-wave phase over the 
cycles in the detector's sensitive bandwidth. The accu- 
mulated phase is 
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where the subscripts in and out denote the values when 
the signal enters and leaves the detector's bandwidth. By 
demanding that the change in &gw caused by the dipole 
term be smaller than tt, one obtains the bound 
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where 77 = fi/m, and M = ri 3 ' 5 m is the "chirp mass" , the 
mass that determines the lowest-order quadrupole effects. 
For LIGO/ VIRGO systems, f m is typically chosen to be 
30 Hz. A more accurate estimate obtained using the for- 
malism of matched filtering (including post-Newtonian 
effects - see below) weakens this bound by about a fac- 
tor of 1.3, assuming a signal-to-noise ratio of 10. The 
resulting bound can be fit by the analytic formula 
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Whether this bound provides a useful constraint on the 
theory depends on the system in question. 

(i) Neutron star and black hole. Since sbh = 0.5, and 
sns < 0.2, S > 0.3 . The resulting bound on cjbd is 
given, from Eq. (|1.5|), by 
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where S/N denotes the signal-to-noise ratio of the de- 
tected signal. The resulting bounds on wbd are plot- 
ted against the black-hole mass, for various neutron-star 
masses, in Fig. [l]. 

(ii) Two neutron stars. For neutron stars, sns varies 
weakly with mass (see for example Table 3 of |1|), so 
that typically S is smaller than 0.05, and for neutron 
stars each around XAMq, is ve ry small indeed. Thus, 
the small value of M. in Eq. ( |l.5| ) is compensated by 
the smallness of S, and the resulting bound on wbd is 
generally weaker than solar-system results unless the dif- 
ference in mass between the two neutron stars exceeds 
about 0.5Mq. For the extreme case of O.7M and 1.4M Q 
neutron stars, the bound could be as large as 1100. The 
inferred bounds are sensitive to the assumed equation of 
state for neutron-star matter. For a particular assump- 
tion about the dependence of S on mass, Fig. || shows 
the bounds that could be achieved, assuming a signal-to- 
noise ratio of 10. 

(iii) Two black holes. Because sbh = 0.5, 5 = 0, 
and there is no dipole radiation at all (see |L5| for dis- 
cussion). In fact the evolution of the system and the 
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resulting gravitational radiation are identical to the gen- 
eral relativistic results, except that the effective gravi- 
tational mass of each black hole is given by Hawking's 
"tensor mass" my, Jl?]] related to the inertial mass by 
rriT = (3 + 2wbd)/(4-I-2wbd)?ti- Since the effective grav- 
itational masses are the only parameters determined from 
the gravitational-wave signal, no test of BD is possible 
from inspiralling two-black-hole systems. 

In order to place a bound on wbd, we must be able 
to decide among cases (i), (ii) and (hi). This requires 
that both the chirp mass M and the reduced mass pa- 
rameter 77 have been measured with sufficient accuracy 
that the mass of one of the bodies is known to be greater 
than the maximum mass for a neutron star, and thus is 
a black hole, and that the mass of the other is known 
to be less than the maximum mass, and is thus likely 
(though not certain) to be a neutron star, or that the 
two masses, if both less than the maximum neutron-star 
mass, are sufficiently different to provide an interesting 
bound. We have extended the matched filtering anal- 
ysis to include post-Newtonian effects in the evolution 
of the gravitational-wave frequency, effects that depend 
explicitly on the reduced mass parameter ?y, and found 
that the accuracy in determining M. and r\ simultane- 
ously with the dipole radiation effect is sufficient for this 
purpose. The dipole radiation effect varies as v~ 2 w r/m 
relative to quadrupole radiation, while post-Newtonian 
correction terms vary as m/r, hence the two effects are 
relatively uncorrelated in the matched filtering (correla- 
tion coefficients of order 0.90). 

The main result is that for a neutron star of mass 
typical of those of well-measured pulsars (1.3 — 1.5M©), 
and a relatively light black hole (3M Q ), a bound on 
the Brans-Dicke parameter around two times the current 
solar-system bound could be obtained. For a low-mass 
system (0.7M Q neutron star, 3M Q black hole), a bound 
around 2000 could be obtained. For double neutron-star 
systems of sufficiently different mass, interesting bounds 
could also result. 

The remainder of this paper provides details. In Sec. 
II, we summarize the relevant BD equations for gravita- 
tional radiation and orbital motion in systems containing 
compact objects. Section III derives bounds on the dipole 
radiation effect using both the crude estimate of accumu- 
lated phase, and a full matched-filtering analysis. In Sec. 
IV we discuss the possibility of observing BD effects in 
the gravitational-wave amplitude, including testing for 
the existence of a third, scalar polarization mode in the 
gravitational waveform. Section V discusses the results. 

II. COMPACT OBJECTS AND GRAVITATIONAL 
RADIATION IN SCALAR-TENSOR GRAVITY 

The equations of motion and of gravitational radiation 
for systems containing compact objects in Brans-Dickc 
theory were first derived by Eardley [fill and extended 



by Will and Zaglauer [ p^JT^ ]. Here we present the basic 
background and formulas needed for our purpose, taken 
from (lf| 

We work to lowest order in an expansion in powers of 
v 2 en m/r, corresponding to Newtonian order for orbital 
motion, and to what in GR would be called quadrupole 
order for gravitational radiation. Yet the equations in- 
clude contributions due to the self-gravitational binding 
energy of the compact objects, contributions that are for- 
mally of post-Newtonian order, but that could be sizable 
because of the relativistic nature of neutron stars or black 
holes. These effects are determined by the "sensitivity" 
of the inertial mass of each body A to changes in the local 
value of the effective gravitational constant G (caused by 
changes in the scalar field). It is defined by 

s A ee -d{\nm A )/d{\nG) . (2.1) 

In the weak-field limit, it is straightforward to show that 
sa is the negative of the usual gravitational self-energy, 
given by sa = (1/2) f A pp'\x — x'| _1 c? 3 a;d 3 a;', where p 
is the mass density. For neutron stars, s depends on 
the equation of state; representative values are given in 
Table 3 of p5| . For instance, for a 1.4M Q neutron star, 
s = 0.143 and 0.125 for equations of state O and M 
respectively. On the other hand, s can vary slowly with 
mass: for equation of state O, for instance, s = 0.1432 for 
to = 1.4OM and s = 0.1440 for to = 1.44M Q . For black 
holes, a dimensional argument shows that m cx G -1 / 2 , 
so sbh = 0.5. 

A. Two-body orbits 

The equation of motion and Kepler's third law for two- 
body orbits with orbital frequency u) (orbital period P) 
take the form 

d 2 x/dt 2 = -Gmx/r 3 , u = 2ir/P = (g m /r 3 ) 1/2 , 

(2.2) 

where 

g = 1 - £(s x + s 2 - 2s x s 2 ) , (2.3a) 

e=(2 + ^ BD ) _1 . (2.3b) 

For a circular orbit, the energy and orbital velocity are 
given by 

E = -(l/2)£/Ltm/r, v 2 = Gm/r . (2.4) 

B. Gravitational radiation energy loss 

The rate of energy loss for a quasi-circular two-body 
orbit (that is, circular, apart from an adiabatic inspiral) 
is given by 
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(2.5) 



where the first term is the combined quadrupole/mono- 
pole contribution, and the second term is the dipole con- 
tribution, and where 



S = si - s 2 , 



r = 1 — 2{m\S2 + m2S\)/m . 



(2.6a) 



(2.6b) 



(2.6c) 



(2.6d) 



The gravitational-waves are dominantly at the frequency 
f = lu/tt, corresponding to twice the orbital frequency. 



C. Gravitational waveforms 



In the far-zone, the spatial components of the radiative 
metric perturbation h^ v = rf v — U (— .g)<? A " / are given by 



(2.7) 



where Greek and Roman indices denote spacetime and 
spatial components, respectively, where ip is the pertur- 
bation of the scalar field <f> about its asymptotic, cosmo- 
logical value <fio and where, to leading order in v 2 ~ m/r, 



9 ij = 2(1 - ^R-^S/df ^mAX 1 ^ 



(4/x/i?)(l - -£)(t/V - Gmx l x J /r 3 ) , (2.8a) 



^M. = C(M/^)|r[(N • v) 2 - gm(N ■ x) 2 /r 3 ] 

-(gT + 2A)m/r-2S(N v)l , (2.8b) 



where R and N are the distance and direction unit vector, 
respectively, of the observer, and A = 1 — si — s 2 + 0(£). 

The components of the Riemann tensor R 0l0 J mea- 
sured by a detector can be shown to be given by R 0l °i = 
— \{d 2 1 ' dt 2 )h % i , where h l i is the effective gravitational 
waveform, given by 



h v 



an 

TT 



1 



(2.9) 



where TT denotes the transverse-traceless projection. 
Note that the full gravitational waveform is transverse 



but not traceless because of the presence of the scalar 
contribution. 

For quasi-circular orbits, the waveform becomes 



R 



Q^ T + S(5 l] - WW) 



fi = 2(l--£) — (\ i X j -n i n j ) 
2 r 



(2.10a) 



(2.10b) 
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4 I 

_25(^) 1 / 2 iV. a|, (2.10c) 
where n = x/r, and A = v/w. 



III. TESTING SCALAR-TENSOR GRAVITY 
USING MATCHED FILTERING OF 
GRAVITATIONAL WAVEFORMS 

A. Phase-shift estimate 

Because broad-band detectors such as the free-mass 
laser interferometric systems detect the gravitational 
waveform h lJ (t) superimposed on the noise, and because 
hundreds to tens of thousands of cycles may be observed 
in the bandwidth, the observations are especially sensi- 
tive to the evolution of th e fre q uenc y and pha se of the 
wave. By combining Eqs. (2J2), ( |2.4| ) and (2J5), one can 
show that the frequency of the waveform evolves accord- 
ing to 



; 96 gVa 



11/2 



5 k d r ^ 2 
96 g m 



(3.1) 



We define the Brans-Dicke chirp mass M. and the dipole 
parameter b according to 



M = ( K 3 / 5 /5 4/ V /5 
b=(5/96)( K - 3/5 g- 6/5 ) KD S 2 . 



(3-2) 



Defining u = nA4f, we put Eq. (3.1) into the form 

u = 7W- 1 (96/5)u 11/3 (l + bT] 2/5 u- 2 / 3 ) . (3.3) 
Integrating, we get 

tt -«/3[i _ (4/5)6t7 2/5 u- 2/3 ] = (256/5)(t c -t)/M, (3.4) 

where t c is the time at which u — > oo . We have expanded 
the expres 
fact that 



the expression (3.4) to first order in brj 2 / b u 2//3 , using the 
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(3.5) 



From Eq. (1.3), the number of cycles observed in a given 
bandwidth can be written $gw = f™° llt (u/u)du, 



giving 



-5/3-1 
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-7/3-, 



(3.6) 



Demanding that the phase contribution of the di- 
pole term be no more than n, we obtain b < 
(1127r/5)?7- 2 / 5 u 1 7 n /3 , where we assume that / out >> f m 
(1000 Hz vs. 10 Hz). To lowest order in 1/ojbb, k = Q = 
I, kd = 2/wbd, and thus b w (5/48)5 2 /wbd, resulting 
in the bound given in Eq. ( |1.4| ). 



B. Matched-filter analysis 

To obtain a more accurate estimate of the bound 
that can be placed on the dipole parameter b, we carry 
out a full matched-filter analysis, following the method 
described by Chernoff and Finn [0] and Cutler and 
Flanagan M. To the accuracy needed, we approximate 
the observed gravitational waveform, Eq. ( 2.10[ ), in a 
given detector by h(t) w 1l{h Q (ty s '^}, where h°(t) is 
the slowly-varying Newtonian-order contribution to the 
waveform amplitude, dependent upon the distance to the 
source, its location on the sky, the orientation of the de- 
tector, and on the source parameters M, rj and r; <fr(t) 
is the gravitational-wave phase, dominantly at twice the 
orbital phase, and TZ denotes the real part. The phase 
includes the dipole and higher order post-Newtonian cor- 
rections, in principle. Calculating the Fourier transform 
of h(t) in the stationary phase approximation, we obtain 



Hf) = { Ar 7/V * 0</</ " 



(3.7) 



where A oc R~ 1 A4 5 / 6 x [function of angles and detec- 
tor orientation], and fmax ~ 0{m~ l ) corresponds to the 
frequency when the inspiral turns into a plunge toward 
coalescence, and 

*(/) = 2vr/i c - $ c - tt/4 

+ ll^ 5/3(1 ~7 6r ' 2/5 ^ 2/3) ' (3 - 8) 
where $ c formally is the gravitational- wave phase at time 

tc- 

With a given noise spectrum S n (f), one defines the 
inner product of signals hi and /12, by 



(hi\h 2 ) = 2 



h\h 2 



Sn(f) 



(3.9) 



The signal-to-noise ratio for a given signal h is given by 

p[h] = S/N[h] = (/i|/i) 1/2 . (3.10) 

If the signal depends on a set of parameters 8 a which are 
to be estimated by the matched filter, then the rms error 
in 8 a in the large S/N limit is given by 



A9 a 



(3.11) 



where S aa is the corresponding component of the inverse 
of the "covariance matrix" or "Fisher information ma- 
trix", T a b, defined by 



ab 



dh 



dh 



(3.12) 



The correlation coefficient between two parameters 9 a 
and 8 b is 



ab 



-^ab , 



'Y]aaJ^bb 



(3.13) 



For a noise spectrum, we adopt the analytic fit to the 
LIGO "advanced detector" noise spectral density, 



Sn(f) 



00 / < 10 Hz, 

So[(/o/7) 4 + 2(l + C/7/o) 2 )] />10Hz, 



(3.14) 



where S =3x 10" 48 Hz _1 and f = 70 Hz. The cutoff 
at 10 Hz corresponds to seismic noise, while the /~ 4 and 
f 2 dependences correspond to thermal and photon shot 
noise respectively. 

We shall adopt the following five parameters to be es- 
timated: In .A, <I> C , fat c , ln.M, and 6, where b = br] 2 / 5 . 
The corresponding partial derivatives of h(f) are 



dhjf) 
d\nA 
dh(f) 
dfotc 
dhjf) 

dhjf) 
dhiM 
dhjf) 
db 



Hf), 

2-Ki{f/h)h{f) , 
-ih(f), 

-^™- 5/3 km-~bu- 2/3 } 

224 w ; 



The signal-to-noise ratio is given by 

p 2 = MA\ 2 f - i/3 I(7)/S , 
where we define the integrals 

- q / 3 { X - 4 + 2 + 2x 2 y 1 dx. 



m = / ■ 

Jl/7 



(3.15a) 
(3.15b) 
(3.15c) 
(3.15d) 
(3.15e) 

(3.16) 

(3.17) 



We also define the coefficients B q = I(q)/I(7). Since our 
a priori expectation is the validity of GR, we are looking 
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for a bound on \b\. Equivalently, we wish to determine the 
error in estimating b about the nominal value b — 0, thus 
we set b = in Eq. ( 3.15d ) . Then the elements of the co- 



variance matrix turn out to be proportional to p 2 u n ^ 3 B q 
for various integers n and q, where uq = irAdfo- In- 
verting the matrix to obtain E ofc , we obtain from the 
elements E MM , £ fcb , and ^ Mb , A(lnM) = 57.6u 5 /3 / p, 

A(b) = 53.2ul /3 /p, and c Mi = -0.98 

Substituting /o = 70 Hz, and working to first order in 
1/wbd, so that M = rf/ 5 m and b = (10/96)7 ? 2 / 5 5 2 /w B d 
(see Eq. (3.2)), we obtain 



A(lnM) = 6.56 x 10" 



M Q J {S/Nj ' 



(3.18a) 



where the final term is the "tail" effect. Adding the pa- 
rameter In rj to the set to be estimated, we find the partial 
derivatives 



dhjf) 
dlnM 
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(3.20b) 
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6.14 x 10- 6 rT 2/5 
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(3.18b) 



C. Post-Newtonian effects 

The interpretation of gravitational-wave observations 
as testing scalar-tensor gravity relies upon determining 
the masses of the two bodies with sufficient accuracy 

(i) that one can decide with confidence whether each 
body is more or less massive than the accepted neutron- 
star maximum mass (modulo the many uncertainties 
in that number) in the neutron-star/black-hole case, or 

(ii) that one can establish that the mass difference ex- 
ceeds some critical value, in the double neutron-star case. 
The determination of the individual masses makes use of 
post-Newtonian corrections to the orbital phase, which 
depend explicitly on the reduced mass parameter r\. It is 
then necessary to check whether such determinations can 
be made simultaneously with the bound on the parameter 
b. To this end, we extend the matched filter analysis to 
include post-Newtonian corrections. Those corrections 
have not been fully calculated to date in the context of 
BD, but one expects them to be the same as those of GR, 
within corrections of 0(1/o>bd)- To incorporate them 
into our analysis in a first approximation, then, it suf- 
fices to add the appropriate GR post-Newtonian terms 
of M, Eq. (3.13), to the Fourier transform phase, Eq. 
(|3.8"|), to obtain 



*(/) = 2nft c - $ c - tt/4 



-5/3 



128 



1 



7 



brj 



2/5 u -2/3 



9 l 336 + 4 m 



T67rf/ 



-3/5, 



(3.19) 



The other partial derivatives in Eqs. (3.15) are un- 
changed. Setting b = 0, we then calculate and invert the 
covariance matrix and evaluate the errors in the five rele- 
vant parameters (the parameter In A decouples from the 
rest, and is not important for our purposes), along with 
the correlation coefficients between A4, r] and b. Notice 
that the fact that b involves 77 implicitly does not affect 
the outcome, because we are considering the results cen- 
tered around b — 0. For various double neutron-star and 
neutron-star/black-hole systems, the results are shown 
in Table J. Note that b is less strongly correlated with 77 
than is A4, a result to be expected because of the very 
different dependences of the dipole and post-Newtonian 
effects on m/r or on u. Nevertheless, a result of the corre- 
lation is to weaken the bound on S 2 /ujbd in Eq. ( 3.18b| ) 
by a factor of about 2.3. Eq. (1.5) provides an analytic 
fit to the result. 

Notice from Table [| that the accuracy in determining 
ij is a few percent, consistent with the results of 0]. To- 
gether with the high accuracy in determining M, this 
will lead to accurate values for the two masses, except 
for the degenerate region around equal masses, rj = 0.25. 
However, that region is excluded from our considerations 
because of the weakness of the bound in the nearly-equal- 
mass, double neutron-star case, and because of the im- 
possibility of determining unambiguously that one of two 
objects of nearly equal mass is a black hole while the 
other is not, in the mixed case. Figure [3] shows the re- 
gions in M. — rj-space corresponding to the three types of 
system. 



D. Dependence on the nature of the system 

The nature of the system being obse rved is important 
to the interpretation of the bound ([h^) as a test of grav- 
itation theory. 

For two neutron stars, the chirp mass is relatively small 
(1.2M Q for two 1.4M Q neutron stars), so the bound on 
iS 2 /wbd can be quite small. The reason is that such sys- 
tems enter the detectors' bandwidth (say, at 10 Hz) with 
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a large separation relative to the total mass (r/m ps 180 
for two 1.4M Q neutron stars). As a result, gravitational 
radiation damping is weaker, and more cycles occur in 
the detectors' bandwidth (up to 16,500 cycles for two 
1.4M Q neutron stars), hence the matched filter is more 
sensitive to effects on the template phase evolution. How- 
ever, the sensitivity difference S is also small, because of 
the weak dependence of the sensitivity on the masses 
of neutron stars, and because the neutron stars are ex- 
pected to have similar masses, as has been seen in known 
binary pulsar systems, such as the Hulse- Taylor system 
PSR 1913+16 Q. For relatively stiff equations of state, 
which are required in order to have neutron stars suffi- 
ciently massive (> 1AM®) to agree with masses inferred 
from known binary pulsars, Table 3 of |lj| shows that 
the sensitivity varies roughly linearly with mass, with 
s ~ O.2771/M0, reaching a maximum near the maximum 
mass for the given equation of state (around 1.46M© for 
the stiff equations of state discussed in Jl5|). To get a 
rough idea of the bounds on wbd that might be possible 
in double neutron-star syste ms, we substitute the relation 
S ps O.2/3<5m/M into Eq. (1.5), where dm = m 1 — m 2 , 



and where f3 may vary between and 1 , and use the fact 
that Sm 2 = (1 — 4i])m 2 , to obtain 



(3.21) 



The resulting bounds for f3 ps 1 are plotted in Fig. [| For 
j3 ps 1, the bounds are weaker than solar-system bounds 
for neutron stars whose masses differ by less than O.5Af . 
At the other extreme, the bound could approach 1000 for 
a 0.7 - 1.4M Q pair. 

For black holes in scalar-tensor gravity, sbh = 0.5. 
This is a consequence of the fact that, in the forma- 
tion of a black hole, the scalar field is radiated away 
(see ]l5| . p^ , |l9l,p0| for further discussion). As a result, 
from Eqs. (O) and (|!|), we see that = 1- £/2, 
k = (1 - £/ 2 rT<S = T = A = 0. Substituting these val- 
ues into Eqs. |0), (EDI), (O) and (O), we see that the 



dipole effects vanish, and the equations become equiv- 
alent to those of GR if we replace every mass by the 
"tensor mass" JItJ , given by 

m T = m(l - £/2) = m(3 + 2cj bd )/(4 + 2w BD ) • (3.22) 

If gravitational radiation is the only information about 
the system, then only the tensor masses are measured 
in a matched filter. Since the behavior of the system is 
general relativistic in terms of tensor masses, then no test 
of BD in double black-hole systems is possible. 

For the case of neutron-star/black-hole systems, 
sbh = 0.5 while Sns ~ 0.1 — 0.2. Although the black 
hole cannot support scalar field in its vicinity, the com- 
panion neutron star can, and consequently dipole gravi- 
tational radiation can occur. For stiff equations of state, 
sns < 0.14, and thus we can approximate S > 0.3. With 
these assumptions, Eq. (1.6) leads to the bounds plotted 
in Fig. | 



IV. TESTING SCALAR-TENSOR GRAVITY 
USING WAVEFORM AMPLITUDE AND 
POLARIZATION 

From Sec. II. C we see that scalar-tensor gravity af- 
fects the amplitude of the gravitational waveform in two 
important ways: (i) it introduces a third, non-traceless 
polarization state (Eq. (2.10a)), and (ii) it introduces, 
through the dipole term, a contribution at the orbital 
frequency in addition to the quadrupole/monopole con- 
tributions at twice the orbital frequency (final term in 
Eq. (2.10c)). However, these are unli kely to be testa ble, 
for the following reason. From Eqs. |[|) and fl3.15| ), it 
is simple to see that the cross-components between In .4 
and the other four parameters in the covariance matrix 
vanish, so that r a fc has the form 



r a h = P 



1 
D 



(4.1) 



where D is a 4 x 4 matrix corresponding to the parame- 
ters $ c , / t c , InM, and b. Thus A(ln^) = (E- 4 - 4 ) 1 / 2 = 
1/p ps 0.1, for a signal-to-noise ratio of 10. How- 
ever, the scalar-tensor corrections to the amplitude are 
much smaller than this. The ratio of the monopole part 
of S ( Eq. J2.10ej )) to the quadrupole term Q« (Eq. 
( gl0t] )) is 0(£/8) fa 2 x 10- 4 (500/w BD ), whil e the ra- 
tio of the dipole part of S to the quadrupole term is 
O^Sir/m) 1 / 2 } « 10^ 3 (500/w BD )(5/0.3)(r/100m) 1 / 2 . 
Thus, despite the presence of qualitatively new contribu- 
tions to the waveform, the lower sensitivity of matched 
filtering to amplitudes, as compared to phases, makes 
those contributions difficult to detect or bound. 



V. DISCUSSION OF RESULTS 

We have shown that interesting, if not spectacular 
bounds on the coupling constant w B d of Brans-Dicke the- 
ory could be obtained by matched filtering of a gravita- 
tional wave signal from an inspiraling binary of a neutron 
star and either a neutron star of very different mass, or 
a low-mass black hole. The bound that can be obtained 
decreases rapidly with decreasing neutron-star mass dif- 
ference and with increasing black-hole mass. 

Another class of laser-interferometric detectors that 
may be relevant to this discussion is the space-based 
class of systems, such as LAGOS @ and SAGITTAR- 
IUS/LISAG (22). These will be most sensitive to gravi- 
tational waves in the frequency range 10 -4 to 10 _1 Hz, 
and could detect gravitational waves from the inspiral of 
a IMq star into a black hole in the 10 6 - 10 7 Mq range. 

Noting that the bound on ujbd scales with the frequency 

-7/3 



/oas/ ~'' , we find from Eq. ( |l.6 ) for a IMq star and a 
10 6 Af Q black hole, with /o ~ 10 -2 Hz, that a bound ex- 
ceeding 50,000 might be possible. This warrants a more 
detailed matched-filter study, including the use of the ap- 
propriate noise curve for space-based detectors analogous 
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to Eq. (3.14), and including the effects of spin of the cen- 
tral hole, which are likely to be important. This study is 
currently under way. 

It is useful to point out that, for generalized scalar- 
tensor (ST) theories, i.e. those in which the coupling 
constant becomes effectively a function of the scalar field 
p3| , the foregoing conclusions still apply, with the fol- 
lowing principal change: the dipole parameter kd, Eq. 
( |2.6b| ), becomes 



K D 



2G 2 



(2 + wst) 



1 + 



2lo!, 



ST 



(3 + 2lo st ) 



(5.1) 



where 



'ST 



dujsT/d(f>. (There are other changes in- 



duced by Wgxi in k j T and other formulae, but they are 
unimportant for our purposes.) In the large-wsT limit, 
this means repla cing wbd by cjst/X1+ w st/2 w |t) 2 m the 
bounds of Eqs. (|L4j), (|l|) and ( pT8b| ). 

Another difference between BD and such generalized 
scalar-tensor theories is that, while in the former case, 
large wbd implies that all physical predictions are close 
to those of GR (roughly within 0(1/wbd)), in the lat- 
ter case the statement holds only for weak scalar-field 
situations (G(f> « 1), such as in the solar system. How- 
ever, in the strong-field interiors of neutron stars, or in 
the early universe, where the scalar field may have values 
very different from its weak-field values, the differences 
between the scalar-tensor theory and GR may be signif- 
icant, despite a large exterior value of c^gx |^]24] , |25| , |2(| . 
In evaluating the sensitivities of neutron stars, we used 
neutron star models computed using GR, since BD is a 
small correction for large wbd throughout the stellar inte- 
rior [^7|. But the strong-field effects in ST theories alter 
the situation, and quite different values of sensitivities 
for a given neutron star mass may result p0| , p6[ . Since 
our quoted bounds were strongly dependent on the sensi- 
tivities, it is impossible to draw strong conclusions about 
bounds on ST theories from inspiralling binaries without 
further research. 

How do the bounds we have suggested compare with 
other empirical bounds and with theoretical expecta- 
tions? Solar-system bounds on scalar-tensor gravity may 
improve from the current level of wbd > 500 during the 
coming decade. The Relativity Gyroscope Experiment 
(Gravity Probe B) aims to measure the geodetic preces- 
sion of a set of orbiting superconducting gyroscopes at 
the 10 -5 level, which would correspond to wbd > 60, 000 
p8| . Orbiting optical interferometers for astrometry may 
measure the deflection of light at the 5 microarcsecond 
level, resulting in bounds on wbd at the level of 10 5 [^9| . 

In BD, cjbd is an arbitrary constant, so no value has 
an a priori theoretical significance (except wbd = oo). 
However, in generalized theories, the value of ujst is cou- 
pled to the dynamics of the scalar field (or fields, in 
multi-scalar theories ^4j), and its present, average cos- 
mic value may depend on the evolution of the universe. 
Damour and Nordtvedt |50| have pointed out that, in 
one class of generalized scalar-tensor theories, in which 



uo{(f)) = — | — i(xln0) _1 , with < (f) < 1, cosmic evolu- 
tion from strongly non-general-relativistic early universes 
tends toward a large-w "attractor" at the present epoch. 
For specific models and values of the parameter \i the 
present value of wst can range from2xl0 4 to 10 6 . Unless 
there are unusual strong-field effects in this theory that 
would modify our conclusions based on BD, these values 
are well above the bounds we anticipate from coalescing 
binaries. 
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FIG. 1. Bounds on ojbd from inspiralling neutron-star/- 
black-hole binaries, plotted against black-hole mass, for var- 
ious neutron-star masses. Hatched portion indicates black 
holes with mass less than 3.OM0, where identification as a 
black-hole may be ambiguous. Curves assume S = 0.3 and a 
signal-to-noise ratio of 10. 



FIG. 2. Bounds on cjbd from inspiralling double neutron- 
-star binaries, plotted against mass of neutron stars. For equal 
masses, dipole radiation is suppressed, and no bound on wbd 
results. Curves assume a linear dependence of sensitivity on 
mass, with S = 0.25m/MQ, and a signal-to-noise ratio of 10. 



FIG. 3. Chirp- mass/reduced-mass parameter plane, show- 
ing location of three different types of compact binaries. 
Chirp mass is plotted in units of the neutron-star maximum 
mass. 
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TABLE I. The rms errors for signal parameters, the corresponding bound on wbd, and the correlation coefficients cmtii c < 



by Eq. (3.14), and a signal-to-noise ratio of 10 is assumed. Masses are in units of M©, At c is in msec. For neutron-star/black 



hole systems, S = 0.3 is assumed; while for double neutron-star systems, S = 0.2Sm/MQ is assumed. 
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